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COMPUTATION  OP  HOLLOW  CYLINDER  EXPLOSIONS 


ABSTRACT 

This  report  treats  the  explosion  caused  by  a  charge  which  Is  placed 
between  two  hollow  co-axial  cylinders.  It  Is  assumed  that  at  the  beginning 
of  the  motion  of  the  cylinders  the  explosive  has  been  detonated  completely 
so  that  the  space  between  the  cylinders  is  filled  by  a  gas  under  high 
pressure.  The  flow  of  this  gas  during  the  collapse  of  the  inner  cylinder 
and  the  expansion  of  the  outer  cylinder  is  governed  by  a  quasi -linear 
second  order  partial  differential  equation.  This  equation  is  solved 
nunerically  using  a  two-level  difference  scheme.  The  solution  furnishes 
the  fragment  velocities  of  the  outer  cylinder  and  the  collapse  velocity 
of  the  inner  cylinder.  The  results  are  compared  with  experimental  data 
and  asymptotic  formulae. 

The  same  computation  scheme  is  applied  to  the  computation  of  the 
motions  of  two  parallel  plates  with  a  charge  exploded  between  them.  The 

results  agree  with  values  obtained  by  approximate  formulae. 

*  .  * 
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INTRODUCTION 


Fragment  velocities  of  exploded  cylinders  have  been  investigated  in 

the  past  mainly  by  experiment.  In  addition  to  experimental  results  some 

approximate  formulae  for  the  velocities  have  been  established  by  Gurney  and 

3  8* 

Sterne  using  energy  considerations.  ’  In  general  these  formulae  provide, 

*  •  •  • 

together  with  adequate  experiments,  satisfactory  results.  However,  if  the 
explosive  is  located  between  two  co-axial  hollow  cylinders,  no  such 
formula  can  be  derived  by  the  method  used  by  Gurney  and  Sterne.  In  such 
cases  one  must  rely  only  on  experiments  which  may  furnish  the  velocities 
of  fragments  of  the  outer  cylinder.  The  measurement  of  the  collapse 
velocity  of  the  inner  cylinder  apparently  requires  more  complicated 
techniques  and  no  publications  of  such  measurements  are  known  to  the 
author.  On  the  other  hand,  the  knowledge  of  the  collapse  velocity  caul 
be  of  some  help  for  studies  of  shaped  charge  problems. 

One  possible  method  to  compute  the  collapse  process  of  the  inner 
cylinder  is  treated  in  this  paper.  We  arrive  at  this  method  by  formulating 
first  an  equation  of  motion  for  the  flow  of  explosive  fumes  in  the  space 
between  the  two  cylinders.  Boundary  conditions  for  this  equation  erne 
established  by  considering  the  properties  of  the  cylinders.  The  solution 

of  this  equation  furnishes  then  the  desired  collapse  and  fragment  velocities. 

%  .  V 

The  equation  of  motion  formulated  emd  treated  here  is  a  quasi -linear 
second  order  partial  differential  equation.  It  is  solved  numerically  by 
using  a  two-level  difference  scheme. 

Supercript  numbers  denote  references  which  may  be  found  on  page  50, 
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Only  some  of  the  results  obtained  by  the  numerical  integration  of  the 
differential  equation  can  be  compared  with  experimental  data,  namely  the 
positions  of  the  fragments  of  the  outer  cylinder  at  given  instants.  This 
comparison  shovs  that  the  computations  are  compatible  with  experiments. 

A  close  agreement  is  obtained  during  the  early  phase  of  the  explosion. 

In  special  cases  the  numerical  solution  can  be  compared  with  the 
approximate  formulae  mentioned  above.  In  all  cases  where  such  formulae 
exist  both  results  agree  completely.  An  equally  close  agreement  of 
results  was  obtained  for  the  case  in  which  the  cylinders  are  degenerated 
into  parallel  plates.  These  agreements  indicate  that  the  numerical 
solution  process  employed  is  sufficiently  accurate. 

Disagreements  between  experimental  results  and  numerical  solutions  may 
be  caused  by  inadequate  assumptions  made  for  the  development  of  the 
theory.  These  assumptions  can  be  changed  and  the  theory  correspondingly 
refined  when  more  accurate  experimental  data  are  available.  The  theory 
developed  permits  such  refinement. 


CHAPTER  1 
PROBLEM  OUTLINE 


Assume  that  the  space  between  two  co-axial  metal  cylinders  is 
filled  with  an  explosive.  The  detonation  of  the  explosive  may  cause 
the  outer  cylinder  to  burst  and  the  inner  cylinder  to  collapse.  We 
are  primarily  interested  in  the  motion  of  the  inner  cylinder  after  such  a 

detonation.  We  obtain  a  similar  problem  If  we  replace  the  cylinders  by 

* 

two  parallel  plates  with  the  explosive  located  between  the  plates. 
Mathematically  this  problem  is  a  limit  case  of  the  first  one  if  we  let  the 
radii  of  the  cylinders  become  infinite. 

The  gas  flow  and  mass  motions  du'ing  the  process  described  depend 
in  a  rather  complicated  manner  on  a  large  number  of  parameters.  In  order 
to  reduce  the  number  of  parameters  we  will  treat  a  simplified  model  of 
the  processes.  This  model  is  obtained  by  assuming  the  following  idealizations 
of  the  real  problem. 

1.  The  cylinders  are  infinitely  long  and  all  motions  are  circular 
symmetric  and  identical  for  all  orthogonal  cross-sections  of 
the  cylinders.  This  assumption  reduces  the  problem  to  that  of 
a  one-dimensional  flow  which  is  cylindrical  for  finite  radii 
of  the  cylinders. 

2.  The  explosive  is  detonated  completely  before  the  cylinders 

•'art  moving.  This  assumption  means  that  the  space  between 
the  cylinders  is  filled  with  a  high  pressure  gas  at  the 
beginning  of  the  motion  of  the  cylinders. 
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rj.  The  gas  between  the  cylinders  is  a  perfect  inviscid  gas 
and  the  gas  flow  is  adiabatic,  i.e. ,  without  heat  exchange 
between  gas  particles.  This  assumption  may  be  Justified  by 
the  observation  that  we  are  interested  in  the  first  phase  of 
the  motion  only,  that  is  until  the  complete  compression  of 
the  liner.  The  flow  during  this  phase  resembles  very  much 
a  simple  wave.  Particularly  we  can  expect  that  no  shocks 
will  develop  in  the  flow  during  this  phase,  the  duration  of 
which  might  be  some  10  seconds. 

4.  The  outer  cylinder  is  surrounded  by  vacuum,  and  vacuum  is 
also  Inside  the  inner  cylinder.  This  assumption  specifies 
the  type  of  the  boundary  conditions  for  the  gas  flow.  A 
complete  formulation  of  the  boundary  conditions  requires 
more  detailed  assumptions  about  the  cylinders.  These 
assumptions  will  be  discussed  in  separate  chapters  treating 
the  boundary  conditions. 

CHAPTER  2 

EQUATION  OF  MOTION 

For  the  description  of  the  problem  we  use  cylindrical  coordinates 
r,  <p,  z  (Figure  2.1).  The  position  of  a  gas  particle  which  was  originally 
at  the  distance  R  from  the  cylinder  axis  is  independent  of  the  coordinates 
cp  and  z  because  of  the  Assumption  1.  Hence  the  status  of  the  system  is 
completely  described  by  the  distance  r  of  any  gas  particle  from  the  axis. 
Our  system  is  thus  one -parametric  and  can  be  described  by  the  function 
r(R,t)  with  the  property 

r(R,0)-R.  (2.1) 
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In  this  Chapter  we  will  derive  a  differential  equation  for  the  function 
r(R,t).  To  this  end  we  consider  first  the  forces  acting  on  a  volume 
element  V  of  the  gas.  At  time  t  =  0  the  element  may  have  the  coordinate  R 
and  its  volume  is  (Figure  2.2) 

V(R,0)  «=  R.dR  hep  +  §(dR)2  hep  (2.2) 


At  time  t  the  gas  which  was  originally  in  V(R,0)  has  moved  to  r(R,t)  and 
it  occupies  the  volume  (Figure  2.3) 


V(R,t)  =  r||  dR-h-ep  +  G(cp.(dR)2) 


(2.3) 


The  forces  acting  on  the  element  are  caused  by  the  pressure  p(R,t) 
within  the  gas.  Because  of  the  circular  symmetry  assumed  we  need  to 
consider  the  resultants  of  the  pressure  in  the  direction  cp  =  0  of  Figure  2.3 
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only.  (Other  components  cancel  each  other.)  These  resultants  are 
F1  =  p.r.h.2.sin  £  cp 

F2  =  F5  =  p||  .dR.h.sin  h  9  +  sin  cp  0((dR)2) 

=  -(p.r  +  p~  dR  +  r|^  dR)  h.2  sin  i  9  +  0((dR)2  sin  £  cp) 


The  sum  of  these  forces  is 

F  =  -  |^.r.dR.h.2  sin  cp  +  0((dR)2  sin  |  cp).  (2.4) 

With  p(R,t)  being  the  density  of  the  gas  at  r(R,t),  we  obtain  from  (2.3) 
and  (2.4) 

+  .  '2-5> 

For  dR  -*  0  and  cp  -«  0  the  last  equation  becomes 


(2.6) 


A  second  equation  we  obtain  by  considering  the  mass  within  the  volume 
element  V.  Since  in  our  system  no  new  mass  is  created,  the  quantity  pV 
must  remain  constant.  We  conclude,  therefore,  frtto  (2.2)  and  (2.3) 

pr||  =  p(R,0).R  +  O(dH) 

or,  with  dR  -»  0 

pr||  =  p(R,0).R.  (2.7) 
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Substitution  of  (2.7)  into  (2.6)  yields 


d^r  1  r  dp 

dt2  =  ■  pfR^y  R*SR  * 


(2.8) 


The  initial  density  distribution  p(R,0)  we  may  consider  as  known. 
However,  the  Equation  (2.8)  still  contains  the  unknown  quantity  dp/dR.  In 
order  to  eliminate  this  quantity  we  have  to  make  some  assumptions  about  the 
properties  of  the  gas.  (So  far  we  have  only  assumed  that  the  gas  is 
lnviscld  in  the  derivation  of  Equation  (2.6)).  Here  we  use  Assumption  5 
which  furnishes  the  following  relation  between  the  pressure  p  and  density  p 
for  an  adiabatic  flow  of  an  lnviscld  perfect  gas: 

(^)  s  0.  (2.9) 

The  exponent  7  is  the  ratio  of  specific  heats  for  the  gas  (Ref.l,p.  10). 
From  (2.9)  we  deduce  the  relation 

p(R,t)  -  p(R,0)«  (2.10) 

Substitution  of  (2.7)  into  (2.10)  furnishes 

P(R,t)  =  p(R,o)»(r  t*ry&)y.  (2.11) 

With  (2.1l)  we  can  eliminate  the  pressure  from  the  Equation  (2.8).  The 
result  is 


at2  =  "  p(R,0)  ’  R  *  3R  {p(R>°)(r<kr/dR)  }  (2.12) 


or 


lU 
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This  final  equation  of  motion  is  a  quasi-linear  second  order  partial 
differential  equation  for  the  function  r(R,t)  in  case  of  a  cylindrical  flow. 
The  functions  p(R#0)  and  p(R,0)  we  consider  as  given  initial  values  for  the 
problem. 

For  comparison  we  derive  the  corresponding  equation  of  motion  for 
cylinders  which  have  degenerated  into  plates.  In  order  to  carry  out  the 
limiting  process  we  write 

r  *=  r  +  R, 

R  =  R  +  R, 

(2.U) 

p(R,t)  =  p(R,t), 
p(R,t)  =  p(R,t) 

and  let  R  -*  «.  The  limiting  process  applied  onto  (2.13)  yields  then  the 
equation 


£*  B  p(R.o)  r ii  .  _ ar.ajli.qn 

p(R,0)  dR  L9r2  p(£,0)  dR  dR  J 


(2.15) 


The  most  significant  difference  between  the  Equations  (2.13)  and 
(2.15)  is  the  factor  r^"^  in  (2.13).  As  the  gas  flow  approaches  the  axis 
of  the  cylinder  this  factor  might  cause  an  increase  of  the  acceleration. 
In  the  equation  of  motion  (2.15)  neither  the  independent  variable  6  nor 
the  dependent  variable  r  appears  explicitly.  This  indicates  the  absence 
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of  any  effects  of  the  geometry  of  the  experiment  onto  the  flow.  For  the 
flow  represented  by  (2.15)  ve  would,  therefore,  expect  a  gradual  decrease 
of  acceleration  to  zero  due  to  the  loss  of  pressure  as  the  gas  expands. 

In  order  to  have  in  the  sequel  a  uniform  formulation  for  cylinders  and 
plates  we  will  use,  whenever  convenient,  for  the  independent  variables 
R  or  l  the  symbol  x  (assuming  a  *  x  £  b)  and  for  the  unknown  functions 
r(R,t)  or  r(R,t)  the  symbol  u(x,t).  Then  the  equation  of  motion  for 
cylindrical  flow  is 


Px(x,°)  i 

VpTxIoT  J 


(2.16) 


and  in  case  of  plates 


u. 


tt 


px(*,°)-i 

ux-p^J* 


(2.17) 


We  see  that  in  both  cases  the  following  three  parameters  determine  the 
flow: 


(1) 

the  ratio  of  specific  heats,  7; 

(2) 

the  function 

-  M 1 

(2.18) 

(3) 

the  function 

P„(x,0) 

®(x)  ‘  "pT^oT  ' 


The  quantity  7  is  a  known  parameter  of  the  explosive.  The  functions  P(x) 
and  Q(x)  depend  in  a  simple  manner  on  the  initial  status  of  the  gas. 
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CHAPTER  3 

SPECIFIC  PROPERTIES  OF  CYLINDRICAL  FLOW 


Investigations  of  centripetal  cylindrical  flows  of  incompressible 

liquid  show  very  high  pressures  and  velocities  when  the  flow  approaches  the 

2 

axis  of  the  cylinder  .  It  has  been  anticipated  that  the  compress¬ 
ibility  of  a  real  liquid  would  reduce  the  pressures  and  velocities 
{loo.  oit.  ).  The  question  arises  whether  such  a  reduction  is  sufficiently 
large  to  cancel  completely  the  increase  of  pressure  in  a  centripetal 
flow.  The  following  investigation  will  demonstrate  that  generally  such 
is  not  the  case  for  adiabatic  flows. 


Assume  that  the  flow  starts  at  time  t  =  0  with  the  following  initial 
values 

r(R,0)  »  R 

"  V<R>  4  °<  (5.1) 

p(R,o)  =  p0(R), 

p(R,o)  *  p0(R). 

R  is,  as  in  Chapter  2,  the  Lagrangian  coordinate  of  the  problem, 
v(R)  is  the  particle  velocity  and  p  and  p  are  the  pressure  and  density, 
respectively. 

By  substitution  of  (3.1)  into  the  equation  of  motion  (2.9)  we  compute 
the  initial  values  of  the  second  time  derivative: 


&  r(R,0) 
dt2 


p;»> 

p^rT 


(5.2) 


For  time  t  >  0  we  obtain  by  expanding  r(R,t),  — and 


in  power  series  of  t 
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(5.3) 


1  2  ,  5, 

r(R,t)  =  R  +  t-v(R)  -  |  td  +  0(tO, 


-  1  +  tv  (R)  -  \  t2<^>  *  0(t}), 


(3-M 


=  v(R)  ' t  p^rT  +  0(t2)* 


(3.5) 


The  corresponding  value  of  p(t,R)  we  obtain  by  substituting  (3.3) 
and  (3.4)  into  (2.1l).  (The  Equation  (2.1l)  was  derived  under  the 
assumption  of  an  adiabatic  flow  from  the  condition  for  conservation  of 


mass.)  The  result  is 


p(R,t)  -  PoGOfrfo  t' 


Po{i  +  t(X-  +  v)+t2(^.ig.i(g)')+tV)}- 
po{1  ■  rt(l  +  v’ )  * 


(5.6) 


♦  ‘W1  <f  ♦  *'>2  -  ^ir  ♦  §  ><r +  £>'>]} +  °<‘5>- 

o  o 

The  difference  between  the  pressure  at  time  t  and  the  initial  pressure 


p(R,t)  -  pq(R)  =  -  P07t(^  +  v’ )  + 


(3.7) 


P*  p' 


+ 1  *  +  x')2  +  -  -')2  +  r +  (r>'] +  0(t5)- 

ro  o 


From  (3*7)  we  conclude  that  the  pressure  in  the  flow  will  increase 
whenever 
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R  v'(R)  <  -  v(R). 


(3.8) 


For  a  centripetal  flow  v(R)  <  0.  Therefore,  (3*8)  can  be  satisfied  for 
any  v'(R)  if  R  is  sufficiently  small.  Hence  we  can  conclude  than  an 
increase  of  pressure  in  a  centripetal  flow  is  possible  even  if  v'  (R)  >  0, 
l.e.,  for  flows  where  the  particles  nearer  to  the  axis  have  higher  velocities. 

In  the  special  case  v(R)  =  const,  the  pressure  will  increase 
throughout  the  flow,  i.e.,  for  all  values  of  R. 


We  consider  now  another  special  case,  namely  pQ(R)  =  const,  and  ask 
whether  the  pressure  gradient  can  become  negative  when  the  pressure  has 
increased.  First  we  note  that  for  Pq(R)  =  const,  the  second  term  in  (3*7) 
is  positive  and,  therefore,  (3.8)  is  sufficient  to  insure  an  increase  of 
pressure  up  to  third  order  terms  in  t. 

The  gradient  of  the  pressure  p(R,t)  is  in  the  case  considered 


dR 


-  P0?t(^  +  v' )'  +  P072t2(^  +  v')(-  +  v' )'  + 
+  P07t2(|  -  v')(|  -  v'  )'  +  0(t5) 


(3.9) 


For  small  t  this  gradient  is  negative  if 

D[v]  =  v"  +  jp  -  \  >  0  (3*10) 

R  R2 

and  any  negative  function  v(R)  which  satisfies  (3*8)  and  (3*10)  will  cause 
an  increase  of  the  pressure  and  a  negative  pressure  gradient.  The 
simplest  example  of  such  a  function  is 

v(r)  •C1  +  CgR.  (3.11) 
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This  function  satisfies  (3*3)  if 

(J  +  v'  )*R  «  Cx  +  2CgR  <  0  (5.12) 

and  It  satisfies  (5.10)  If 

(|  -  v’)R  *  cx  <  0.  (5.13) 

Clearly,  if  (3*12)  and  (3*13)  axe  both  satisfied  the  function  v(R)  <  0. 
Furthermore,  this  function  makes  the  second  term  in  (3*9)  negative  and 
the  third  term  equal  to  zero.  Therefore,  positive  terms  in  the  expression 
(5*9)  for  the  gradient  of  the  pressure  are  at  least  of  third  order  in  t. 


Figure  5.1 

For  comparison  we  carry  out  the  same  analysis  as  above  for  cylinders 
degenerated  into  plates.  Because  the  condition  for  the  conservation  of 
mass  is  now  different  we  obtain  instead  of  (3.6)  the  following  equation 


for  the  pressure: 
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P  I 

=  p0(R){l  +  tv'  (R)  -  |  t2(^)’  +  0(t3)} 
a  P0(R){l  -  7  t  V'  (R)  + 


(3.14) 


p 1 

+  \  +  i)v'2  -  ^(p2)’]}  +  0(t5). 


As  a  condition  for  the  increase  of  pressure  we  obtain  instead  of  (5*8) 
simply 

v' (R)  <  0.  (3.15) 

The  pressure  gradient  of  an  initially  constant  pressure  becomes  now 

"  -  r  t  ▼"  +  +  i)v'v"  +  o(t5) .  (3.16) 

A  negative  gradient  of  the  increased  pressure  will  therefore  be  obtained 
if  v(R)  satisfies  in  addition  to  (3*15)  the  condition 

v"(R)  >  0.  (3.17) 


As  in  the  case  of  a  cylindrical  flow  both  conditions  together  insure  the 
increase  of  pressure  and  negative  gradient  up  to  terms  of  third  order  in  t. 

Figure  3.2  shows  a  typical  velocity  profile  satisfying  (3*15)  and 
(3.17).  The  main  differences  between  the  profiles  shown  in  Figures  3*1 

* 

v(r)  satisfies  in 
(0,R)  the  conditions 
v'  <  0  and  v"  >  0. 

Figure  3*2 
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and  3*2  is  that  the  pressure  in  a  cylindrical  centripetal  flow  increases 
ndt  only  in  all  such  cases  where  it  would  increase  in  a  linear  flow  but 
also  in  other  cases  with  positive  velocity  gradient.  (Condition  (3*15) 
implies  (3*8)  because  v  <  0.)  Therefore,  pressure  increase  in  a 
cylindrical  centripetal  adiabatic  flow  is  possible  and,  if  the  velocity 
gradient  remains  finite,  likely  for  small  R.  The  maximum  pressures  obtained 
in  such  flows  depend  on  the  boundary  conditions  which  have  not  been 
discussed  in  this  analysis. 

CHAPTER  4 

BOUNDARY  CONDITIONS  FOR  THE  CASING 
Boundary  conditions  for  the  gas  flow  at  the  outer  cylinder  can  be 
established  by  considering  the  forces  acting  on  a  cylinder  element  shown 
in  figure  4.1.  Thereby  we  need  to  consider  force  components  in  the 
direction  cp  =  0  only,  because  the  other  components  cancel  each  other  due  to 


The  resultant  of  the  pressure  in  the  direction  cp  =  0  is 

F1  '  P(R2't)’r(R2,t)’h’2  8in  \  (^-l) 

where  h  is  the  height  of  the  cylinder  element. 
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About  the  stresses  within  the  metal  we  assume  that  they  can  be 
approximated  by  an  average  value  a  across  the  thickness  Tg.  Assuming  a 
linear  elasticity  law  with  the  modulus  of  elasticity  E  we  obtain  then 
the  following  expression  for  the  force  component  Fg  resulting  from  stresses 


F 


2 


o  Tgh  2  sin  ^  cp  = 

(4.2) 

r(R  ,t)  -  Rp  , 

E - - - -  -Tg-h-2  sin  |  cp. 


The  mass  of  the  cylinder  element  is 

M  =  p2  R2  Tg  h  cp  +  p2  i  Tg  h  cp 


(4.3) 


where  p2  is  the  density  of  the  cylinder  material. 

Equilibrium  of  forces  acting  on  the  cylinder  element  furnish  the 
equation 


p2  B2  T2(1  +  £jl*S 


.  2  sin  ^  cp 

-E-5TT2>  - 


r-R, 

C 

*2 


(4.4) 


For  cp  -*  0  we  obtain 


d2r  r(R2,t)  p  (R2,t)  E  r(R2,t)  -  Rg 

^2  =  p2  T2  R2(l  +  T2/2R2)  '  ‘  R2(l  +  T2/2R2)*  (4*5) 

The  pressure  p(Rg,t)  can  be  expressed  in  an  adiabatic  flow  with  (2.11)  by 

R  7 

p(R2,t)  =  p(R2>°)(r(R2#t).ar(R2,t)/dR)  •  (4.6) 

Substitution  of  (4.6)  into  (4.5)  yields  then  the  boundary  condition  for 
R  =  Rg.  Expressed  in  terms  of  the  variables  x  and  u(x,t)  this  boundary 


condition  is 
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m 

4 


« w(b,t)  =  p<b)gft*2l .  ft7'1 1_  ■  .  a 

tt'  ’  P«W0  'u'  y  p~b 


2  2 


U*  P2' 
X 


T 

-  b  2 
b  *  W2> 


(4.7) 


where  P(x)  =  p(x,0)/p(x,0)  is  a  parameter  of  the  differential  equation 
(see  (2.18))  and  Wg  is  defined  by 


W2*T2<1+^)- 


(4.8) 


The  quantity  W2  has  the  dimension  of  a  length  and  it  characterizes 
together  with  p2  the  casing.  Normally  T2  <  <  R2  and  consequently, 

Wg  w  T2,  i.e.,  Wg  is  approximately  equal  to  the  thickness  of  the  casing. 

The  boun^  ry  condition  (4.7)  is  of  the  form 


utt(b,t)  -  Fm  -  Fg, 


(4.9) 


i.e.,  the  acceleration  of  the  casing  is  proportional  to  the  difference 
between  a  term  which  depends  on  the  mass  of  the  casing  and  a  term 
Fg  which  depends  on  the  strength  of  the  casing. 


Fg  has  the  form  indicated  in  (4.7)  during  the  first  phase  of  the 


movement  only.  Once  the  yield  stress  a  in  the  cylinder  is  surpassed, 

y 

the  linear  stress-strain  law  ought  to  be  replaced  by  a  different  assumption. 
Assuming  for  instance  a  constant  stress  in  the  cylinder  for  the  phase 
between  the  yield  point  and  rupture  point  we  obtain  for  Fg  the  approximation 


/  _E_  u(b,t)  -  b  _2 
P2b  ‘  b 


w 


u  <  b(l  +  o  /E) 

y 


F 


S 


=  < 


V 


if  b(l  +  ay/E)  2  u  <  b(l  +  er)  (4.10) 

if  b(l  +  ej  2  u  . 

r 
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wr 


The  quantity  e  in  (4.10)  is  the  deformation  of  the  cylinder  at  the 

3  8 

rupture  point.  Experiments  have  indicated  that  for  steel  cylinders  * 

e  0.2  -  0.6.  The  formula  (4.10)  for  F0  may  be  replaced  by  a  more 
r  o 

complicated  formula  if  a  contraction  of  the  casing  during  the  movement 
is  possible.  Generally  in  such  cases  a  hysteresis  in  the  stress-strain 
law  must  be  taken  into  account.  The  inclusion  of  such  conditions  in 
(4.7)  is  straightforward  and  most  easily  done  at  the  level  of  numerical 
calculations . 

At  the  beginning  of  the  motion  Fg  =  0  and  FM  is  large.  As  u(x2,t) 

increases,  Fg  increases  too,  whereas  F^  decreases.  In  all  technically 

interesting  cases,  however,  F^  is  at  the  rupture  point  still  larger  than 

F0  by  several  orders.  Therefore,  the  term  F-  might  be  dropped  in  (4.7) 
o  o 

altogether.  On  the  other  hand,  its  inclusion  into  the  boundary  conditions 

is  straightforward  and  may  be  done  anyhow,  especially  if  the  computations 

are  done  by  fixed  machine  codes.  Its  inclusion  in  such  codes  will  save 

the  user  from  the  necessity  to  check  the  significance  of  Fg  in  any  particular 

case.  The  parameters  E,  0  and  c  of  F_  are  normally  known  more  accurate 

y  r  o 

than  the  parameter  P(b)  in  the  term  FM» 

The  first  term  F^  in  the  boundary  condition  (4.9)  represents  the 
action  of  the  pressure  of  gas  on  the  inner  surface  of  the  casing.  After 
rupture  of  the  cylinder  the  gas  may  escape  freely  through  the  gaps  and  it 
has  been  suggested  that  no  further  acceleration  of  the  cylinder  fragments 
takes  placed  The  latter  assumption  will  probably  underestimate 
the  fragment  velocity.  In  order  to  obtain  an  upper  bound  for  that  velocity, 
we  may  assume  that  the  fragments  retain  their  shape  after  fragmentation 

and  are  accelerated  by  a  force  proportional  to  the  gas  pressure  times  the 
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surface  area  of  the  fragments.  With  this  assumption  we  obtain  for  the 
formula 


/.  .  7-1  ,  x  (1  +  e  ) 

P(b)  ^  <*)  i  0>-U> 


M 


In  case  the  cylinders  are  degenerated  into  plates,  no  considerations 
about  the  strength  of  the  material  are  necessary.  The  right  side  of  the 
boundary  condition  (4.9)  is  then  given  simply  by 


Fm  =  P(b) 
Fs  =° 


Pfao) 

P2T2 


ux(b,t)-7 


(4.12) 


Note,  that  in  this  case  the  product  PgTg  enters  the  boundary  conditions 
instead  of  P2W2’ 

CliAPTER  5 

BOUNDARY  CONDITIONS  FOR  THE  LINER 

Boundary  conditions  for  the  gas  flow  at  the  inner  cylinder  ("liner") 
can  be  established  by  considerations  similar  to  those  of  Chapter  4.  Let 
the  density  of  the  cylinder  material  be  and  the  initial  thickness  of  the 
cylinder  be  T^.  The  outer  initial  radius  of  the  liner  is  R^.  As  boundary 
conditions  during  a  first  phase  of  the  motion  we  obtain  then  an  equation 
of  the  same  type  as  (4.9) >  namely 

<Wa't>  =  fM  -  V  (5-1) 

The  first  term  on  the  right  hand  side  of  (5*1)  represents  the  action  of  the 
pressure  and  it  depends  on  the  mass  of  the  cylinder: 
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with 


(5-5) 


w^  is  a  parameter  of  the  liner,  approximately  equal  to  itr  thickness. 
P(x)  is  the  parameter  of  the  differential  equation  defined  by  (2.18). 


The  more  interesting  and  problematic  second  term  in  (5*1) ,  which 
depends  on  the  strength  of  the  material  will  be  discussed  further  below. 


If  the  explosion  is  sufficiently  strong  to  cause  the  liner  to  collapse, 
then  eventually  it  will  be  compressed  into  a  solid  cylinder  which  we  assume 
is  rigid.  Consequently  the  boundary  condition  changes  into 


u(a,t)  =  umin  =  const.  (5*^) 

after  this  stage  is  reached.  The  value  of  u.  is 


u 


min 


1  max 


2  w1R1  , 


(5-5) 


where  p,  is  the  density  of  the  liner  after  total  compression. 

1  max 

Normally  for  metal  cylinders  one  would  assume  p.^  max  =  p^. 

The  complete  boundary  conditions  are  then 

utt(a,t)  ■  *„  -  *g  »<*»*>  >  W  (5>6) 

ut(a,t)  =  0  for  u(a,t)  =  u^. 

The  term  fg  in  (5*6)  can  be  computed  for  the  beginning  of  the  motion 
by  assuming  a  linear  stress-strain  law.  We  obtain  thus  for  the  elastic 
phase 
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(5.7) 


u(a,t)  -  a 
a 


a(l 


o/E)  <  u(a,t)  5 

v 


a 


where  E  is  the  modulus  of  elasticity  of  the  liner  and  is  its  yield 


stress. 


Little  is  known  about  the  behavior  of  the  liner  if  after  the  yield 
point  is  surpassed  further  compression  takes  place.  Obviously  the  liner 
will  be  buckled  or  even  damaged  during  this  phase.  Fragments  of  material 
will  be  driven  toward  the  axis  of  the  cylinder  and  pushed  into  a  wedge. 
They  may  be  pressed  together  into  a  new  layer  immediately  after 
fragmentation  adding  thus  to  the  strength  of  the  cylinder.  It  is  as  easy 
to  give  plausible  reasons  for  an  increase  of  fg  as  for  a  decrease  of  fg 
during  this  phase.  As  a  working  hypothesis  we  may,  therefore,  assume  that 
fg  ■  const,  after  the  yield  point  is  surpassed.  The  complete  formula  for 
fg  is  then 


_E_  .  u(a,t)  -..a.J:  if  a(l  _  0^E)  <  u(a#t)  s  a 

fg  ='  W  (5*^ 

JL  .  JL  .  -2  if  Vn  <  S  a(l  "  av/E) 

p^a  E  w1  y 

As  for  the  casing,  we  can  easily  include  a  hysteresis  type  stress-strain 
relation  in  the  formula  for  fg  if  expansion  of  the  liner  during  the 
motion  takes  place. 


The  hypothesis  about  fg  can  be  checked  by  measuring  the  liner 
velocity  in  experiments.  The  function  fg(u(a,t))  can  then  be  modified 
to  fit  the  experimental  results.  An  indirect  testing  of  the  working 
hypothesis  is  possible,  too.  We  may  assume  that  the  work  required  for  the 
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deformation  of  the  liner  is  transformed  into  heat.  The  compressed  liner 
should  have,  therefore,  a  high  temperature  after  the  explosion.  (Heating 

of  the  liner  by  heat  conduction  should  be  negligible  because  of  the  short 

-5 

duration  of  the  process.  The  compression  is  completed  within  some  10 
seconds^.)  The  heat  contents  of  the  liner  can  be  measured  and 
compared  with  the  theoretical  value  according  to  the  working  hypothesis. 
With  the  assumption  (5.8)  the  work  against  the  stresses  becomes 

W  =  TT  h  Tx  Rx  oy  {i  +  j~  (uyleld  -  unln)},  (5.9 

where  u  is  the  radius  of  the  liner  at  the  yield  point  and  h  is  the 

yield 

height  of  the  cylinder.  The  first  term  in  the  braces  in  (5-9)  represents 
the  contribution  to  the  work  during  the  elastic  phase  of  compression  and 
the  second  term  is  the  contribution  due  to  the  working  hypothesis. 

The  temperature  increase  of  the  liner  due  to  W  is 


T  = 


W 


Mlcm 


(5.10) 


where  is  the  mass  of  the  liner  and  ca  is  its  specific  heat. 
Substituting  (5.9)  into  (5.10)  we  obtain  for  the  temperature  increase 


ff. 


(5.H) 


Ir.  case  the  liner  is  degenerated  into  a  plate  the  boundary  condition 
becomes  particularly  simple.  In  this  case  we  have 


fM  ■  -  P(a)p^rf  ‘  ux(a't)'T  • 

fs*°- 


(5.12) 
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CHAPTER  6 
INITIAL  CONDITIONS 


We  obtain  simple  initial  conditions  for  the  problem  if  we  assume  that 
a  constant  volume  detonation  has  turned  the  explosive  into  gas  completely 
some  time  before  the  cylinders  start  moving.  We  can  then  assume  that  the 
pressure  and  density  are  constant  throughout  the  gas  and  that  the  gas  is 
at  rest.  The  parameters  of  the  differential  equation  (see  (2.18)) 
become 


P(lt)  =  3  Po  *  con8t-’ 


«(x) 


Pv(x,0) 

“pTxToT  s  °' 


(6.1) 


and  the  initial  values  for  u(x,t)  are 

u(x,0)  =  x 
ut(x,0)  h  o 


(6.2) 


The  numerical  value  of  Pq  is  easily  obtained  from  the  energy  Eq  per 
mass  unit  of  the  explosive.  After  the  constant  volume  detonation 
the  whole  energy  released  is  transformed  into  heat.  Since  we  have  assumed 
that  the  explosive  fumes  are  a  perfect  gas  we  obtain  the  equation 
(Ref.  l,p.  9) 


or 


E 


o 


_  JL_  p(x,0)  _  1  p 

7-1  p(x,0)  ~  7  -  1  o 


(6.5) 


P  =  (7-1)  •  E  . 
0  w  '  o 


(6.4) 


The  values  of  Eq  and  7  are  known  quantities  fo?  a  given  explosive. 
For  example,  for  pentolite  we  have  the  values 
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6  2  2 

E  =  1.152  [kcal/g]  =  4.822 *10°  [m  /sec  ], 

7  =  3.2, 

P  =  1.06.10^  [m2/sec2]. 
o 

If  the  space  between  the  cylinders  is  filled  only  partly  with 
explosive,  the  initial  value  p(x,0)  will  be  a  fraction  of  the  explosive's 
density.  However,  the  pressure  p(x,0)  is  reduced  correspondingly  in  such 
cases  and  the  parameter  Pq  remains  unchanged.  Hence  this  parameter  is  a 
characteristic  property  of  the  explosive,  independent  of  the  particular 
loading,  as  demonstrated  by  (6.4). 

The  loading  density  p(x,0)  enters  the  boundary  conditions  only,  as 
shown  by  (4.11)  and  (5*2). 

The  assumption  of  a  constant  pressure  and  density  throughout  the  gas 
at  the  beginning  of  the  process  is  certainly  a  quite  rough  approximation 
to  the  reality.  On  the  other  hand  it  is  virtually  impossible  to  furnish 
such  non-constant  generally  valid  initial  conditions  which  may  be  justified 
by  the  setup  of  experiments.  The  initiation  of  the  detonation  is  in  a  real 
case  invariably  done  at  one  or  both  ends  of  the  cylinder  and,  therefore, 
the  real  problem  is  not  one -dimensional  anyhow.  A  constant  initial  status 
seems  not  to  be  a  more  serious  idealization  than  the  other  assumptions 
outlined  in  Chapter  1.  The  final  Justification  of  the  idealizations  is, 
of  course,  the  experiment.  If  the  numerical  results  obtained  by  solving 
the  differential  Equations  (2.16)  or  (2. IT)  do  not  agree  within  reasonable 
errors  with  experiments,  we  have  to  reconsider  all  idealizations. 
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CHAPTER  T 

NUMERICAL  SOLUTION 


The  technical  problem  outlined  in  Chapter  1  was  formulated  in 
Chapter  2  as  a  quasi  •linear  second  order  partial  differential  equation. 
The  initial  and  boundary  conditions  for  the  equation  were  derived  in 
Chapters  4,  5  end  6.  For  the  numerical  solution  of  this  equation  we  will 
propose  in  this  Chapter  a  two-level  difference  scheme.  Such  a  scheme  uses 
data  at  t  *  t^  only  to  compute  the  solution  at  t^+1  =  t^  +  1^^.  Usually, 
for  second  order  problems  difference  schemes  with  more  levels  are  used. 
The  reasons  for  applying  a  two-level  scheme  on  the  present  problem  are 
as  follows: 


(1)  A  two-level  scheme  permits  a  simple  change  of  the  size  of 
the  time  step  1^+1,  depending  on  computed  results.  This 
property  is  of  great  practical  value  if  dealing  with  non-linear 
equations. 

(2)  The  solution  by  the  two-level  scheme  presented  in  this 
Chapter  furnishes  concurrently  with  u(x,t)  also  ut(x,t). 

Since  for  the  present  problem  the  values  of  at 
boundaries  are  of  principal  interest,  the  scheme  meets 
practical  needs  by  avoiding  a  numerical  differentiation  of 
the  results. 


The  difference  scheme  will  be  described  only  briefly  in  this  Chapter. 

A  detailed  analysis  of  the  scheme  will  be  published  in  a  forthcoming  paper. 

The  difference  scheme  is  applied  to  a  rectangular  grid  in  the 
x,t -plane  defined  by 
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=  x^  ih  ,  (i  =  0,1,2, . . .  ,n) , 
k 

85  tQ  +  £  ^  >  (k  =  °>1>2f‘*10  =  °)' 

J=0 

(see  Figure  7.1).  The  correct  values  of  u(x,t)  at  the  grid  points  we 
denote  by 


Figure  7.1* 

In  order  to  reduce  the  number  of  indices,  we  will  denote  in  the  sequel 
the  partial  derivatives  d/dt  by  dots  and  the  partial  derivatives  d/dx  by 
apostrophes . 


The  basis  for  the  difference  scheme  are  the  relations 


Ui,k+1  =  ui,k  +  *k+l  ui,k  +  5  ^+1  Ui,k  +  Z  \+l  ui,k+l  +  R1  ^7*2^ 


with 
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(7.3) 


vith 


.  ljj 

R1  *  "  u  ~W>  (0  *  ®1  S  1} 


Ui,k+1  “  Ui,k  +  2  ^k+l  Ui,k  +  2  \+l  ui,k+l  +  R2 


R2  *  -  uU*(x1,tk-»e2lk+1)  (OS02S1).  (7-5) 


For  the  approximate  solution  we  drop  the  remainder  terms  and  R2  in  (7*2) 
and  (7.4),  respectively.  In  order  to  compute  U  and  U  at  grid  points  with 
t  *  t^+^  by  the  truncated  formulae  we  need  then  the  values  of  U,  U  and  U 
at  t  «  t^  and  t)  at  t  =  t^+^.  The  first  two  quantities  are  known  for  the 
first  time  step  from  initial  conditions.  Proceeding  with  the  computations 
hy  one  time  step  at  a  time  we  have  these  values  available  for  t  =  tfc  when 
cooputation  of  U  and  0  for  t  =  t^+^,  starts.  The  value  of  U  we  obtain  by 
using  the  differential  equation 


U  =  d(x,t,U,U,U')  U”  +  e(x,t,U,U,U' ). 


(7.6) 


The  values  of  the  arguments  x,  t,  U  and  U  of  d  and  e  are  known  for 
t  «=  t^.  The  values  of  U'  and  UM  are  computed  by  numerical  differentiation 
of  U.  The  necessary  formulae  follow  from  the  relations 


with 


ui,k  "  2h  (ui+l,k  “  ui-l,k^  +  S1 


\  -  -  uw(Xi  +  e5t,tk)|-  ,  (-1  £  e5  S  1) 


(7.7) 


(7.8) 


ui,k  =  2  (ui+l,k  "  ^i^k  +  ui-l,k^  +  S2 


(7-9) 


with 


S2  -  -  uIV(xi  +  6,h,tk)  •  Ig  ,  (-1  ^  \  ^  (T.10) 

For  the  numerical  differentiation  of  U  we  drop  S1  and  S£  in  (7-8)  and 
(7.10),  respectively. 

The  last  quantity  needed  in  (7*2)  and  (7*M  is  ^ 

approximated  as  follows: 


First  we  compute 


1  ,2 


i,k+l  =  Ui,k  +  Ui,k  +  2  \+l  Ui,k, 
i,k+l  =  \k  +  *k+l  %k' 


*U 


*Ui,k+l  “  2h  ^#Ui+l,k+l  "  *Ui-l,k+l^ 


*uI,k+i  =  h  <*ui+i,  k+i  -  ^ui,k+i +  *ui-i,k+1}- 

'  h 


(7.H) 

(7.12) 

(7.15) 

(7.H0 


Using  these  values  and  the  differential  Equation  (7*6)  with  t  ®  \+i»  we 
obtain  then  a  value  *Uljk+1*  This  value  we  use  as  approximation  of 

Ui,k+1  in  (7*2)  '7*U)* 

By  this  all  terms  on  the  right  hand  of  (7*2)  and  (7*4)  are  estimated 

(except  for  the  remainder  terms  Rx  and  R£  which  are  dropped) .  The 

trunctated  formulae  (7-2)  and  (7-U)  thus  furnish  the  values  of  ^i>k+1  and 

U  and  the  computing  process  starts  anew  with  k  +  1  =  k. 

1  ,k+l 

It  can  be  shown  that  the  approximations  obtained  by  this  computing 
scheme  have  the  following  error  orders  if  u^k  and  u^k  are  known  exactly 
and  certain  higher  order  derivative f  of  u  are  bounded: 
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Error  order  of  Ulk+1:  0(j£+1)  +  0(l|+1h'  ), 

Error  order  of  Ul  k+1:  0(l£+1)  +  oO^h8). 

We  note  that  these  error  orders  are  the  same  if  the  computational 
scheme  is  applied  to  the  more  general  non-linear  differential 
equation  u  *  f(x,t,u,u,u*  ,u") . 

The  flexibility  of  the  scheme  with  respect  to  changes  of  the  time 
interval  1  cm  be  exploited  by  choosing  at  every  time  step  a  different 
interval  in  accordance  with  some  criterion.  A  reasonable  choice  from 
the  view  point  of  the  physical  phenomena  involved  is  to  select  1^+1  such 
that  the  domain  of  dependence  for  the  point  x^,  tk+1  is  for  t  =  t^  within 
the  interval  x1_1  *  x  *  x1+1  (Ref.  5»  P-  275)*  Because  the  domain  of 
dependence  can  depend  on  x  we  chose  for  1^+1  the  corresponding  minimum 
over  all  x^: 

W-f?  (T.15) 

• • • ,n-i 

An  investigation  of  the  stability  of  the  scheme  (in  the  sense  of 
Richtnyer^)  shows  that  the  amplification  matrix  satisfies  the  necessary 
von  Neumann  stability  condition  if 


W 


(7.16) 


Hence,  by  chosing  1^+1  in  accordance  with  (7*15)  we  remain  well  apart 


from  the  von  Neumann  instability  bound. 
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CHAPTER  8 

TESTS  OF  THE  SOLUTION 

In  order  to  test  the  feasibility  of  the  idealizations  assumed  and 
of  the  numerical  solution  process  used,  some  gas  flows  were  computed  and 
compared  with  experiments  and  asymptotic  formulae. 

The  only  experimental  data  available  were  contained  in  a  report 
about  fragment  velocities  obtained  from  positions  of  fragments  of  the 
casing  at  different  time  intervals  after  explosion  .  Unfortunately  this 
report  does  not  give  the  particular  positions  observed  but  gives  certain 
averages  and  average  velocities  instead.  Therefore,  these  data  can  be 
used  to  compare  the  orders  and  trends  of  the  fragment  velocities  only. 

A  description  of  the  cylinders  exploded  in  the  experiments  is  given 
in  Table  I.  The  fragment  positions  for  Experiments  2  and  3  and  the 
corresponding  results  of  computations  are  shown  in  Figure  8.1. 

The  computed  values  u(b,t)  are  in  both  examples  plotted  in  Figure  8.1 
higher  than  those  of  the  experiments,  except  fo*  the  first  parts  of  the 
curves.  Such  a  general  behavior  is  to  be  exp  •  jted  because  the  computation 
does  not  take  into  account  the  finite  lengths  of  the  cylinders  which 
permit  the  gas  to  escape  through  the  ends  of  the  cylinders.  Moreover,  the 
fragment  positions  plotted  are  obtained  by  averaging  the  positions  of 
fragments  from  ends  of  the  cylinder  and  from  the  central  zone  of  the  cylinder. 
For  the  present  purpose  a  comparison  with  the  higher  values  from  the  central 
zone  would  be  more  reasonable.  A  third  effect  of  unknown  magnitude  is  the 
effect  of  the  surrounding  air  on  the  fragment  positions. 
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Table  I.  Parameters  of  Experiments 


Height 
of  Cylinder 
h 

CM  us  C\l 

NO  CM  UN 

t—  un  on 
•  •  • 
4  Oi  O 

Pi 
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Density 
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CO  cO  oo 
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Eh  H  fa 


CM 
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O  to 

«  S 

to  o 
to  .0 

i  o 

•H 
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The  report  about  the  experiments  does  not  give  fragment  positions  for 
Experiment  1  where  two  hollow  cylinders  were  involved.  Averages  of 
the  ratios  v  =  (distance  flown  by  fragments) /(flight  time)  for  various 
z o  les  are  reported  instead.  These  values  are  compared  with  the  computations 
in  Table  II.  The  larger  discrepancy  between  observed  and  computed  values 
for  long  flight  times  is  certainly  to  a  great  extent  due  to  the  fact  that 
the  exploded  cylinder's  diameter  is  almost  twice  as  large  as  its  height. 
Therefore,  a  close  agreement  with  the  behavior  of  an  infinite  long 
cylinder  cannot  be  expected.  Note  also  the  large  variation  of  average 
v  -  values  in  the  experiment. 

In  all  three  experiments  a  closer  agreement  with  computed  values  is 
obtained  for  short  times  after  the  explosion.  The  agreement  becomes  less 
good  with  increased  time  as  unknown  factors  influence  the  fragment 
velocities.  Since  the  main  objective  of  the  computations  is  the  study  of 
the  behavior  of  the  cylinders  during  very  short  times  after  the  explosion 
(i.e.,  during  the  collapse  of  the  liner),  the  comparison  with  experiments 
can  be  considered  as  satisfactory. 

For  Experiments  2  and  3  the  fragment  velocities  can  be  computed  also 
by  approximate  formulae  derived  by  Gurney  and  Sterne  and  quoted  in  the 
appendix.  These  formulae  are  derived  under  the  assumptions  that  the 
inner  cylinder  is  rigid,  the  density  of  the  gas  is  constant  and  the 
gradient  of  the  particle  velocity  of  the  gas  is  constant,  too.  With  these 
assumptions  the  velocity  of  the  outer  cylinder  (i.e.,  fragments)  can  be 
computed  equalizing  the  kinetical  energy  of  the  system  with  the  energy 
contained  in  the  explosive.  Since  the  assumptions  of  constant  density 


1*0 


Table  II.  Comparison  of  Experiment  1  with  Computations 


and  constant  velocity  gradient  are  certainly  most  closely  satisfied  after 
considerable  expansion  of  the  cylinders,  the  velocities  computed  with 
Sterne's  formulae  can  be  considered  as  approximations  to  the  asymptotic 
values  of  the  fragment  velocities.  A  comparison  of  Sterne's  values  with 
those  obtained  by  the  solution  of  the  differential  equation  can  be 
considered,  therefore,  as  a  check  of  the  numerical  integration  process. 

The  values  obtained  are  compared  in  Table  III.  The  close  agreement  between 
these  values  excludes  the  numerical  integration  as  a  reason  for  differences 
between  experiments  and  calculations  for  large  t.  These  differences  are, 
therefore,  caused  by  the  idealizations  of  the  problem. 

Some  typical  results  of  the  calculations  are  shown  in  Figures  8.2,  8.3 
and  8.4.  The  Figure  8.2  shows  some  pressure  profiles  of  the  gas  flow 
corresponding  to  Experiment  1  at  various  instants.  It  is  interesting  to 
note  that  the  pressure  increases  within  the  centripetal  flow  near  the  axis 
already  efore  the  liner  is  compressed  into  a  solid  cylinder.  Such  a 
behavior  of  the  flow  was  anticipated  in  Chapter  3.  After  the  compression 
of  the  liner  the  flow  is  reflected  at  the  now  rigid  inner  cylinder  and  a 
high  pressure  wave  proceeds  into  the  gas.  The  pressure  decays  rapidly 
with  increasing  radius  because  of  the  geometry  of  the  flow.  At  the  time 
the  wave  reaches  the  casing  the  pressure  has  dropped  below  2000  atmospheres 
in  this  example. 

The  Figure  8.3,  tipper  part,  shows  the  positions  of  the  liner  and 
of  the  casing  of  Experiment  1  as  functions  of  time.  It  permits,  however, 
the  following  different  interpretation.  Assume  that  the  detonation  was 
initiated  at  the  right  hand  end  of  a  cylinder  and  that  it  proceeds  with 
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Table  III.  Comparison  of  Asymptotic  Formulae  with  Results 

of  Integration 


Fragment  Velocities  [m/sec] 

Experiment 

Number 

Asymptotic  Value 

Result  of  Integration 

2 

1865  for  Urupt  =  1.6-Rg 
1900  for  urupt  =  1.8-Ra 

1878  at  t  =  113  V.  sec 

3 

1852 

1845  at  t  =  170  m-  sec 

Sandwich 

VA  '  1535 

vB  -  3308 

V.  =  1336 

at  t  =  104  [i  sec 

vB  =  3308 

The  thicknesses  and  densities  of  the  plate*  in  the  "Sandwich"  case 
were  assumed  equal  to  the  thicknesses  and  densities  of  the  cylinders 
in  Experiment  1. 


Figure  8.2  Pressure  profiles 


pttonation  front 
(V  «  7  800  m/«ic ) 


Figure  8.3.  Collapse  profile 
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a  constant  speed  to  the  left.  We  may  consider  the  Figiire  8.3  then  as  a 
picture  of  the  partly  damaged  cylinder  at  a  fixed  time.  The  lower  part 
of  the  Figure  shows  such  an  interpretation  with  the  assumption  that  the 
vf1.)city  of  the  detonation  front  is  78OO  m/sec.  It  would  be  highly 
i  n  'cresting  and  instructive  to  compare  the  computed  collapse  profile  of 
the  liner  with  experimental  data. 

Figure  8.4  is  identical  to  the  upper  part  of  Figure  8.3  with  the 
paths  of  gas  particles  included.  The  form  of  these  paths  suggests  the 
presence  of  a  shock  caused  by  reflection  of  the  gas  flow  at  the  rigidly 
compressed  liner.  Since  this  shock  develops  only  after  the  liner  is 
completely  compressed  it  does  not  influence  the  behavior  of  the  liner 
during  collapse.  From  a  practical  view  point  there  is,  therefore,  no 
need  to  employ  a  more  sophisticated  computing  process  which  permits  to 
consider  entropy  changes  caused  by  the  shock. 

In  the  present  case,  however,  the  computations  were  carried  on  after 
the  development  of  the  shock  in  order  to  obtain  data  which  could  be 
compared  with  experiments.  A  discussion  of  the  effects  of  the  neglected 
entropy  changes  on  these  data  is  therefore  in  place. 

The  most  striking  behavior  of  the  high  pressure  wave  in  Figure  8.1 
is  its  rapid  decay  which  is  caused  by  the  geometry  of  the  flow.  At  the 
time  the  wave  reaches  the  casing  the  pressure  has  dropped  to  a  6mall 
fraction  of  the  initial  pressure.  Since  the  geometry  affects  a  shock 
in  a  similar  manner  we  can  expect  that  an  outward  moving  shock  will 
also  rapidly  decay.  A  small  shock  can  be  replaced  for  computational 
purposes  by  a  wave\  Therefore,  the  difference  between  shock  and  wave 


calculations  becomes  smaller  with  Increasing  time  and  distance  from  the 
axis.  Tb£  magnitude  of  the  pressure  In  the  wave  is  such  that  Its 
Influence  on  the  casing  Is  negligible  in  the  example  computed.  In  this 
example  which  is  typical  for  the  applications  we  have  moreover  other 
more  serious  effects  on  the  behavior  of  the  casing  neglected,  namely  the 
effect  of  the  fragmentation.  The  fragmentation  of  the  casing  takes 
place  somewhere  between  the  points  A  and  B  in  Figure  8.4.  The  shock 
(or  wave)  reaches  the  casing  after  or  Immediately  before  fragmentation. 

The  position  of  the  casing  before  fragmentation  is,  therefore,  not 
Influenced  by  the  reflected  shock.  Any  comparison  with  experiments  should 
be  made  during  this  phase.  Once  the  casing  is  fragmented  the  agreement 
between  computed  and  measured  positions  is  likely  to  became  bad,  whatever 
the  accuracy  of  the  computation  of  the  flow.  The  computation  of  a  shock 
would  influence  the  results  (position  of  the  casing)  only  for  this  phase. 
Hence  even  for  large  time  values  (as  in  the  present  case)  the  possible 
merits  of  a  consideration  of  the  shock  are  questionable. 

In  order  to  test  the  computations  for  cylinders  which  are 
degenerated  into  plates  the  computational  results  were  compared  with 
Sterne's  formulae  for  sandwich  flows.  Experimental  data  for  such  flows 
were  not  available.  For  the  plates  the  thicknesses  and  densities  of 
the  cylinders  of  Experiment  1  were  assumed  and  also  the  same  lor  ling 
density  vas  used.  The  resulting  velocities  are  listed  in  Table  III  and 
they  agree  perfectly  with  Sterne's  values. 
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CHAPTER  9 
CONCLUSIONS 

The  computation  of  hollow  cylinder  explosions  by  the  process 
outlined  in  this  paper  furnish  reasonable  results  in  spite  of  the  coarse 
idealizations  assumed.  The  fair  agreement  of  computed  results  with 
experiments  indicates  that  a  refined  theory  lacks  a  testing  basis  at 
the  present  time.  Therefore,  a  more  sophisticated  theoretical  treatment 
of  the  problem  may  be  delayed  until  more  precise  and  more  detailed 
experimental  results  are  made  available.  In  Chapters  5  and  8  some  special 
experiments  are  proposed  for  this  purpose. 

The  method  presented  provides  simple  means  to  compute  fragment  and 
collapse  velocities  for  long  co-axial  cylinders.  The  accuracy  of  the 
computations  cannot  be  checked  thoroughly  because  corresponding 
experimental,  data  are  not  available.  Available  data  agree  with 
computations  for  times  of  the  order  of  10  tisec  after  the  explosion. 
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APPENDIX 


STERNE'S  FORMULAE 

In  Chapter  8  approximate  formulae  for  the  fragment  velocities  were 

used.  These  formulae  are  derived  by  Sterne^.  A  survey  about  other 

0 

formulae  of  this  type  and  their  derivation  is  given  by  Henry  .  Here  we 
quote  only  the  two  Sterne '  s  formulae  used . 

Case  1.  Hollow  Steel  Cylinder  with  Rigid  Core. 

Let 


Eq  ■  energy  per  unit  mass  of  the  explosive, 

A  =  ratio  of  outer  cylinder  mass  by  explosive's  mass, 


R2  =  initial  radius  of  the  outer  cylinder. 


R  =  radius  of  the  outer  cylinder  of  rupture.  Sterne^  gives 
rup  3 

the  value  R  =  1.6'Rg  whereas  Henry  gives  the 


value  R 


rupture 


1.2-R2, 


R  *  radius  of  the  rigid  core, 
core 


The  fragment  velocity  is  then  given  by 


A  + 


2E 


1* 
Z~ R 


o _ 

rup 

rup 


+  R 
+  R 


core 

core 


For  Rcore  *  0  we  obtain  the  corresponding  formula  for  a  cylinder 
filled  solidly  with  explosive. 

Case  2.  Explosive  between  two  Parallel  Plates.  (Sandwich  Flow.) 
Let 


Eq  =  energy  per  unit  mass  of  the  explosive, 

A  =  ratio  of  the  mass  of  plate  "A"  by  explosive's  mass, 
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V 


B  =  ratio  of  the  mass  of  plate  "B"  by  explosive's  mass 

v  v  =  respective  velocities  of  the  plates  "A"  and  "B". 

A*  B 

The  velocities  are  given  by 

v l  -  (wA+BHi^B+iaEy (1+2B)2’ 

V®  -  (l+A+B)(2jA+l.l)+iaE7  (1+2A)2- 

I 

I 


1 
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